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ABSTRACT 

Bekenstein's (2004) TeVeS theory has added an interesting twist to the search for dark 
matter and dark energy, modifying the landscape of gravity-related astronomy day by 
day. Built bottom-up rather than top-down as most gravity theories, TeVeS-like the- 
ories are healthily rooted on empirical facts, hence immediately passing sanity checks 
on galaxy rotation curves, solar system constraints, even bullet cluster of galaxies 
and cosmology with the help of 2eV neutrinos. Nonetheless, empirical checks are far 
from perfect and complete, and groups of different expertises are rapidly increasing 
the number of falsifiable properties of the theory. The theory has also been made 
much simpler and more general thanks to the work of Zlosnik, Ferreira, Starkman 
(astro-ph/0606039, 0607411). Here I attempt a tutorial of how to compute lensing 
convergence, time delays etc in TeVeS-like theories for non-spherical lenses. I gave 
examples to illustrate a few common caveats of Dark-Matter- guided intuitions. 
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1 INTRODUCTION 



1.1 GL in any co- variant metric theory 

While gravitational lensing (GL) is a cornerstone of Ein- 
steinian gravity, light propagation is actually well-defined in 
any generic metric theories of gravity and one can test such 
theories using lensing data. In fact, light bending is a gen- 
eral property of E&M waves propagating following Fermat's 
principle, which happens in a non-uniform medium where 
the effective speed of light c varies (e.g., as in atomospher- 
ical seeing) in a flat space-time. Lensing also happens as 
light with constant speed c following geodicics in vaccum in 
a curved space-time (e.g., bending by the Sun), but the link 
between the matter density and the curvature need not be 
given Einstein's equation. 

For example, the gravitational pre-factor G e g needs not 
be a true constant in galaxies where the gravity is so weak 



The solar gravity on Pluto 4 x 10 4 greater than a typical 
place in a galaxy. E.g., the Sun's acceleration around the 
Galaxy 



(200 kms" 1 ) 2 
lOkpc 



0.1c 



lm lAngstrom 



HubbleTime day 2 



(1) 



A mundane example of lAngstrom per second squared 
gravity is the mutual Newtonian gravity of two nearly paral- 
lel sheets of printing papers approximately. It is well-known 
that a parallel plate electric capacitor yields a different E- 
field if immersed in vaccum entirely (Casmir effect), or filled 
in a dielectric air inside. The gravitational attraction of two 



sheets of paper could depend on enviornment in a similar 
way for very different physics. Consider a Gedanken experi- 
ment with a gravitationally torquing pendulum made by two 
misaligned suspended sheets of paper. If one could measure 
the period of the torquing pendumlum not only here on 
Earth (as in free-fall experiments in an Einstein tower), but 
also take the table-top experiments to the edge of the solar 
system (where Pioneer 10/11 probes are), in the interstellar 
space (where galactic stars orbit) and in the expanding void 
between galaxies, then one could measure how G e ff changes 
with space and time. 



1.2 The three pillars of the standard cosmology 

The standard cosmological paradigm is built on three pillars: 
Cold Dark Matter, a cosmological constant, and Einsteinian 
gravity. While independent experimental basis of each of the 
three is debatable on astronomical scales, but their synergy 
(characterised by the cosmological pie) has proven amaz- 
ingly successful at describing the Universe especially on large 
scale. 

Despite its apparantly enticing simplicity, the paradigm 
has much to be understood and is under pressure to be 
modified by observations of galaxy scale and by competing 
theories like TeVeS. For example, the experimentally un- 
detected cold dark matter (generally thought to be Super- 
Symmetry particles) is predicted to clump in scale-free fash- 
ion, while observations of dwarf galaxies suggest a kpc-scale 
free-streaming length of dark matter particles. The idea of 
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introducing a constant or a vaccum energy 100 orders of 
magnitude lower than what the SUSY physics can provide 
naturally is still regarded by many theoreticians as unsatis- 
factory. 

These fine-tunings have lead some to believe the 
paradigm is an effective theory, e.g., a 4D projection of a 
more fundamental 5D brane world theory. Some also ques- 
tion the Einsteinian gravity since its associated equivalence 
principles, remain untested on galaxy scale and cosmological 
scale. 

1.3 Modified gravity: motivations and history 

Modifying gravity is a reoccuring excercise which started 
ever since the general acceptance of Einsteinian gravity, 
which was itself a revolutionary modification to Newtonian 
gravity. Many theories modify the Einstein-Hilbert action 
to introduce a new scalar field which manifests itself only 
through the extra bending of space time, but its coupling to 
the metric is different from the simple coupling of massive 
particles with the space-time metric. 

By construction the theories would respect Special Rel- 
ativity prescription of metric co- variance, and preserve con- 
servations of momentum and energy. They do allow for 
a table-top Cavendish-type experiment with a torquing 
pendulum to measure an effective gravitational constant 
G e ff(t,x) which varies with time and enviornment of the 
experiment. For example, the recent motivation to replace 
the cosmological constant in General Relativity leads to the- 
ories with G e ff depending on the curvature of space-time, 
which evolves with the cosmic time in a way to drive the 
acceleration of the universe at late time. 

However, among two dozen theories proposed after GR, 
very few survive the precise tests on SEP in the solar system 
and the well-studied binary pulars. Even fewer are motivated 
and succeeded in addressing both astronomical dark matter 
and cosmological constant. 

1.4 The secret of TeVeS and variants 

TeVeS is an exception. It holds the promise of explaining 
both dark matter and cosmological constant by relaxing the 
SEP (strong equivalence principle) only in untested weak 
gravity envionments like in galaxies, but respecting the SEP 
to high accuracy in the solar system. 

Crudely speaking such theory has an aether-like field 
with an aquadratic kinetic term in its Lagrangian density, 
so the G e f f can be made a function of the strength of gravity 
\g\, such that G e ff is constant within 10 -16 anywhere in the 
solar system, yet varies by a factor of 10 in galaxies and in 
the universe over a Hubble time where \g\ is much smaller. 
Enhancing the G e f / mimics the effects of adding dark mat- 
ter, and reducing the G e ff can drive the acceleration of the 
universe. 

2 A CHARACTERISTIC SCALE FOR DARK 
MATTER 

As one of the important issue to be understood about dark 
matter, it has long been noted that on galaxy scales dark 



matter and baryonic matter (stars plus gas) have a remark- 
able correlation, and respects a mysterious acceleration scale 
a (Milgrom 1983, McGaugh 2005). 

The Newtonian gravity of the known matter (baryons 
etc.) gfc and the dark matter gravity guAf are correlated 
through an empirical relation (Zhao and Famaey 2006, An- 
gus, Famaey, Zhao 2006, Famaey, Gianfranco, Bruneton, 
Zhao 2006) such that the light-to-dark ratio, experimentally 
determined to fit rotation curves, satisfies 

fffc = gpM + agk ^ 

9dm 0,0 

where < a < 1 is a parameter. Let a = we get a very 
simple relation 

gDM ~ y/g~kao, ao = lAngstromsec~ 2 (3) 

where ao is the forementioned gravity scale, below which 
DM and DE phenomena start to surface. 

Such a tight correlation is difficult to understand in a 
galaxy formation theory where dark matter and baryons 
interactions enjoy a huge degrees of freedom. This spiral 
galaxy based empirical relation is also consistent with some 
elliptical galaxies and gravitational lenses. 



3 A SCALE FOR DARK ENERGY 

Equally peculiar is the amplitude of vaccum energy density 
A, which is of order 10 120 times smaller than its natural 
scale. It is hard to explain from fundamental physics why 
vaccum energy starts to dominate the Universe density only 
at the present epoch, hence marking the present as the turn- 
ing point for the universe from de-acceleration to accelera- 
tion. 

This is related to the fact that 

a ~ \/A ~ cH (4) 

where ao is the characteristic scale of DM as well. 

Somehow dark energy and dark matter are tuned to shift 

a 2 

dominance when the energy density falls below g^y- These 
empirical facts should not be completely treated as random 
coincidences of the fundamental parameters of the universe. 
The explanation with standard paradigm has been unsatis- 
factory. 



4 THE METRIC AND DYNAMICS OF THE 
TEVES FIELDS 

TeVeS, as GR, is a metric theory. Let </ M „ being the physical 
coordinates, then near a quasi-static system like a galaxy, 
the physical space-time is only slightly curved, and can be 
written as in a rectangular coordinate x — (xi,xi,X3) cen- 
tred on the galaxy as 

—c 2 dr 2 = gttdt 2 + g rr dl 2 , (5) 
dl 2 — (dx 2 + dx 2 + dx%). (6) 

Introducing a small quantity ^1 < 1, we can write the 
metric components 

9rr « -c g tt w 1 - -J". (7) 
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To show that takes the meaning of a gravitational 

potential, we note that a non-relativistic massive particle 
moving in this metric follows the geodesic 



d 2 Xi 
IF 



dgu 
2dxi 



0. 



d 2 x . , 

^-«-V*(x), 



(8) 



which is the equation of motion in the non-relativistic limit 
where 



4.1 Vector or scalar 

Near a quasi-static system like a galaxy, goo = — (1 + 2<J>), 
where we omit the factor c 2 for clarity. TeVeS predicts a 
time-like vector field with four components, which are ap- 
proximated as 



A a = (1 - 0-$, 0,00) 
and 

A a = -(!-<£ + $,0,00) 



(9) 



(10) 



to the lowest order, where <j> is a scalar field. 

For most of the system that we are interested, the key 
is the module of A, which is equivalently described by the 
scalar field <j) related through 



A 2 = g ab A a A b = -e 



-20 



< 



(11) 



This shows the vector field is more fundamental than 
the scalar field in TeVeS and TeVeS is can be described the 
physical metric and vector field alone. The original proposal 
of Bekenstein contains two metric, while most recent work of 
Zlonik, Ferreira, Starkman (2006, PRD. 74, 0404037) shows 
that the theory is equally described by a single physical met- 
ric g M „, whose geodicics particles and light will follow. The 
other metric (called Einstein metric) is fully described once 
the vector field is specified. 



4.2 TeVeS as dark matter 

In TeVeS, the galaxy potential $ comes from two parts, 
$ = $ fe „ + <^ (12) 

where the known Newtonian gravitational potential <E>fc(x) 
of known matter of density pk (x) satisfies 



V • V$ fcn = 4irGp kn 

and the added scalar field satisfies 

V [p s \7(f>] = 47rGp fc „, 
where 



Lis = 



ao 



C)( 



0,0 



(13) 



(14) 



(15) 



The picture to keep in mind is that the scalar field re- 
places the usual role of the potential of the Dark Matter. 
The vector field A is fully specified once 4> and $ are given. 

To illustrate how the scalar equation come from La- 
grangian of the vector field theory, let's consider a toy model. 



5 AN E&M-LIKE 4- VECTOR POTENTIAL 

Zlonik, Ferreira, Starkman (2006, astro-ph/0607411) gener- 
alised TeVeS as part of a broader class of Einstein- Aether 
theories, which we will follow here. We will neglect a certain 
Lagrangian multipler for normalisation of vector field and 
also make simplifications where possible (setting ci = C3 = 
— 1 and C2 = in their notation, and choosing the simplest 
modification to the Langrangian) . We emphasize the simi- 
larity of the 4-vector field here with the 4-vector potential 
(Ao, A\, A2, A3) in electromagnetism. 

Let the vector field A a — gapA 13 be a time-like unit 
vector with the constrain equation 



A a A 



-1 



(16) 



and let it be coupled to the metric g a g, so that the system 
is governed by an action S or Lagrangian density L, 



-I 



j4 

a x\ 



R 



16ttG 



+ L kn + 1 - 



f + o(f: 



where R the usual Ricci scalar of the metric g, Lkn is the 
known matter Lagrangian density coupling the known mat- 
ter with metric but not coupled to the vector field directly. 
The 3rd term is the local Lagrangian density of the vector 
field, which apart from a dielectric-like modification factor 
(1 — 2//3 + 0(/ 2 )), is the normal electromagnetism-like ki- 
netic coupling to metric 



4f 



32-kG ~ 32ttG ^ 
where the dimensionless / parameter is a measure of the 
strength of the Maxwell tensor field F a/ 3: 



d a A g - d g A 



(19) 



which is the covariant derivative of the 4-potential A a , (as 
for the electric and magnetic field in electromagnetism). 



5.1 Einstein eq. and vector field eq. of motion 

Taking variations of the action respect to the metric g a (i 
and A a respectively we get the gravitational field equations 
and vector equation of motion for this theory respectively. 
Combining the two we have 

Gg/3 

8^G 



Tin, lL I r~f~\ 
a/3 + J-af3, 



(20) 
(21) 



where the left-hand side is the proportional to the Einstein 
tensor Gnu = R^v — ^g^u and on rhs the 1st term is the 
stress-energy tensor of known matter, the 1st term is the 
stress-energy tensor for the vector field T a p, 

T a p = f a0 + ^A a VMf + O(f))g^F„ ], (22) 

which is a non- linear function of derivatives of the field Apy 
Note how the vector field creates the effect of additional mat- 
ter. 

Near a galaxy, Goo = 2VV$, so the 00th moment of 
the above equation reduces to 



V(^i s V$) = 4irGp kn , 
where 



(23) 
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f + 0(f), f: 



fflo 



(24) 



This way we recover the classical MOND equation of Beken- 
stein and Milgrom (1984) in the weak field limit f — > 0, i.e., 
the gravity V3> drops as s/GMao/r far away from a point 
mass M. 



5.2 Hubble expansion equation 

For homogeneous flat cosmology, we can set the metric 
ds 2 = -c 2 dt 2 + a{tf{dr 2 + r 2 d0 2 + r 2 sin 6 2 d(f> 2 ) (25) 



and the vector field 
A a = -A a = (1,0,0,0), 
and its derivatives 



(26) 



/ = 0. (27) 

The Einstein equation reduces to the following eq. for Hub- 
ble expansion 



/ da \ 2 _ 8ttG 



\adtj 



Pkn- 



(28) 



Interestingly this is not modified from the GR, where p is 
the density of known matter of the uniform background. 

The actual theory of TeVeS and Zlonik et al. is more 
sophisticated than above toy model. E.g., Zhao & Famaey 
(2006) proposed to modify TeVeS with a p, = //(l + /), 

2,2 2 

hence the Lagrangian ^(l-^d^g = (/-ln(l + /)) T ^ = 
(l-^+0(/ 2 ))L/, such that we recover Einstein-Newtonian 
in strong gravity, and achieve a good fit to rotation curves 
in the intermediate regime, and have a smooth transition 
between galaxies and cosmology. The term f 2 or Lf can 
contain derivatives other than the Maxwell tensors. It is also 
possible to modify the homogeneous cosmology in TeVeS- 
like theory to create the effect of dark energy. The cosmology 
prescribes distance-redshift relations, essential for lensing. 
For these we assume the background cosmology in TeVeS- 
like theories must fit the same SNe distance data, hence the 
resulting distance-redshift relation (at least at redshifts of 
order unity) is exactly as in LCDM. A detailed study to this 
effect is shown in Zhao (2006, astro-ph/0610925). 



6 THE INTERPOLATING FUNCTIONS IN 
MOND VS TEVES 

In the classical MOND a spherical galaxy of rotation curves 
V(r) could be explained by modifying gravity: 



9 = 



V 2 



GM 



r 2 p,(g/a ) ' 



(29) 



where p,(g/ao) is the effective "dielectric constant", depend- 
ing on the gravitational field strength, g. The standard in- 
terpolating function 



p(x) 



(30) 



is often used in fitting rotation curves. But Zhao & Famaey 
(2006) argued that this function has undesirable features in 
TeVeS. Instead they (Angus, Famaey, Zhao 2006) proposed 
to use 



H(x) 



2x 



1 + (2 - ax) + - ax) 2 + 4:X 



(31) 



which is a parametric a-family, which recovers Bekenstein's 

(2004) toy model and the simple model of Famaey & Binney 

(2005) by setting a = and a = 1. 

The gravitational potential in MOND theory satisfies a 
modified Poisson's equation, 



V[/iV$] = V 2 $jv = ^Gp kn 



(32) 



where the ptn is the density of all known matter. This is 
different from TeVeS, where the total potential is the sum 
of Newtonian potential ($jv) and a potential due to a scalar 
field (fa): 



(33) 



We can see that the scalar potential plays the role of the dark 
matter gravitational potential, and the Poisson-like equa- 
tions for the scalar field relates it to the Newtonian potential 
<I>jv (generated by the baryonic matter), 



V[/i 3 V0 3 ] = V 2 $iv = 4ttGp 



(34) 



where jj, s is a function of the scalar field strength g s = | V0 S |, 
and is derived from a free function in the action of the scalar 
field. In spherical symmetry, we have 



p s g s = n(g s +5iv) = 9n 



(35) 



where /i is the interpolating function of MOND, hence the 
two interpolation function is related by 



(Is 



/' 



(36) 



7 LIGHT BENDING IN SLIGHTLY CURVED 
SPACE TIME 

Light rays trace the null geodesies of the space time met- 
ric. Lensing, or the trajectories of light rays in general, are 
uniquely specified once the metric is given. In this sense light 
bending works ^exactly* the same way in any relativistic 
theory as in GR. 

Near a quasi-static system like a galaxy, the physical 
space-time is only slightly curved. Consider lensing by the 
galactic potential $(r). A light ray moving with a constant 
speed c inside follows the null geodesies dt — yj—^j^dl. An 
observed light ray travels a proper distance l os — l ts + l ot 
from a source to the lens and then to an observer. Hence 
it arrives after a time interval (seen by an observer at rest 
with respect to the lens) 



dt 



dl 

c 



2$(r) dl 



(37) 



containing a geometric term and a Shapiro time delay term 
due to the <E> potential of a galaxy. 

In fact, gravitational lensing in TeVeS recovers many 
familiar results of Einstein gravity in (non-)spherical geome- 
tries. Especially an observer at redshift z — sees a delay 
At obs in the light arrival time due to a thin deflector at 
z = Zi 



cAt obs (R) 
(1 + Zi) ~ 2D t Di 



(R-R s f 



dl 



2$(R,/) 



(38) 
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as in GR for a weak-field thin lens, $/c 2 <§C 1. A light 
ray penetrates the lens with a nearly straight line seg- 
ment (within the thickness of the lens) with the 2-D coor- 
dinate, R — D18, perpendicular to the sky, where Di(zi) = 
'oi/(l + Zl) is the angular diameter distance of the lens at 
redshift zi, D a is the angular distances to the source, and 
Di a is the angular distance from the lens to the source. The 
usual lens equation can be obtained from the gradient of the 
arrival time surface with respect to R. i.e., 

DiDi s 



D s 



-a x (x,y) = x 3 



(39) 



DiDis , , 

y — y> — a v( x 'V) = y° 

where the deflection 

,2d x $(x,y,l) 



f 

a x = I 

J —i. 

f 

J — < 



tu- 



rn 



di 



2d y $(x,y,l) 



and the convergence 



2D 3 



(d x a x + d y a y ) . 



and likewise for the shear, 

72 = D t d y a x , 

and 

7i = -y (d x a x - d y a y ) 
and for the amplification A 
A' 1 = (1 - k) 2 - 71 - 72. 



(41) 



(42) 



(43) 



(44) 



The time delay between a pair of images i and j is given by 
the path integral 



dpL{p), 



(45) 



L(P) = K ^+^ -cW) 

+ ^f»«(P*) + »«fe)_ ay(p) 
dp V 2 

where x(p),y(p) defines a path from image i to image j as 
p varies from pi to . 



8 SUBTLE DIFFERENCES OF TEVES 
LENSING AND DM LENSING 

8.1 k of TeVeS and DM 

An interesting point is that in GR k is proprotional to the 
projected surface density of known matter. This is not the 
case for a non-linear theory. We can express k into of the 
critical density as follows, 



t{x,y) 



AirGDiDis 



(47) 



where we define an effective projected density as follows, 



E(x,2/): 



dlp(x,y,l), 



(48) 



note the integrand is NOT the true matter volume density 
at (x,y,l), rather 



P(x,y,l) 



v 2 Hx,y,i) 



Pkn + peDM > Pkn 



(49) 



4ttG 

because $ is addition of two fields, and we have an effective 
Dark Matter (eDM) from the field, 

V 2 <f)(x,y,l) 



PeDM = 



4ttG 



(50) 



The eDM tracks the known matter, because the TeVeS <j> 
field is determined by non- linearly with pk- 

_ VfV4>(x,y,l) _ V 2 [$-0Qr,jy,Q] 

E.g., a TeVeS point lens has a non-zero convergence due 
to the non-zero effective DM halo. So important differences 
between lensing in TeVeS and in GR is in the predicted 
metric or potential $ for a given galaxy mass distribution 
pk- 

In general, our non-linear Poisson equation can be 
solved by adapting the code of Ciotti, Nipoti, Pasquale 
(2006). In some cases, one can also take the potential-to- 
dcnsity, and start with a reasonable guess for the potential, 
and find the density by taking appropriate derivatives, e.g., 
the application of Angus et al. (2006a, 2006b) on the bul- 
let cluster. In special cases, one can also solve the TeVeS 
Poisson equations analytically. This is the case for Kuzmin 
disks. 



8.2 k in a non-spherical model 

Here we illustrate how to model non-spherical lens galaxy in 
TeVeS, and point out some subtle difference with DM. We 
present models with identical rotation curves, one in TeVeS 
and one in DM, demonstrate the subtle difference in their 
lensing signal. 

To keep the calculations tractable, we approximate the 
lens potential as that of an edge-on Kuzmin disk with a 
density 



p(R,Z) = 



MbS(Z) 



2n(R 2 + ft 2 ) 3 / 2 



(52) 



in cylindrical coordinates around the symmetry axis (R, Z) 
of the Kuzmin disk of total mass M. We work out the bend- 
ing angle and the lens equation in the sky directions x and 
y. Here the Newtonian gravity g is given by, 



JV JVn 

Jr,9z) 



GM 



[R 2 + {Z + 6±) 2 ] 3/2 



(R,Z + b ± ) 



(53) 



where b± = ±6 for positive or negative Z respectively, be- 
cause the upper part of a Kuzmin disk is in a spherical po- 
tential like that of a point-mass centered on (xo,yo, «o) be- 
low the disk. The potential of the lower part of the Kuzmin 
disk is a mirr or-image of th e upper part. The Kuzmin scale- 
length b — J x 2 + y 2 + 2 ( 2 is smaller than the half-light ra- 
dius Thai f by a factor v3- 

The nice thing about a Kuzmin disk is that the 
MONDian gravity is parallel to the Newtonian gravity 
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riigourously because the equal potential contours coin- 
cide with equal Newtonian gravity contours. Assuming the 
Bekenstein mu-function, or 



ao 



we get the TeVeS scalar field 

-(9r,9z) 
where 



r 2 + (z + b±y 



:(R,Z + b±) 



2 
V 



VGMa . 
One can easily verify that 



g 



4irGp, 



(54) 



(55) 



(56) 



(57) 



One can also replace the scalar field with a spherical 
DM halo potential «g/21n(7? 2 + Z 2 + b 2 ) with a gravity 



<9 



DM DM\ 
R ,9z ) 



R 2 + Z 2 + b 2 



(R,Z) 



(58) 



such that they two descriptions give identical gravity ga on 
the equtorial Z = plane where rotation curves are mea- 
sured. So the circular velocity curve is 



V c l = Rg R (R,0) == 



GMR 



[R 2 + b 2 } 



3/2 



VpR 2 
R 2 +b 2 



(59) 



which is rising as solid-body at centre and asymptotically 
flat velocity vo at large radii. Far away from the disk plane, 
TeVeS resembles a spherical DM halo very well (Read & 
Moore 2005). 

Nonetheless there is a subtle difference. While the DM 
force is centrally-pointing, the scalar force is pointed more 
to a point below/above the disk. At large radii 7?, and small 
\Z\, we have 



1.9 



DM I 

\gz I 



b\/GMa Q 



IZI < b < R 



(60) 



R 2 + b 2 ' 

so the TeVeS scalar field provides *stronger vertical force* 
than the spherical DM close to the disk plane. 



8.3 Lensing by Kuzmin Disk 

Assuming Bekenstein's /x function, we can get the lens equa- 
tions, 



_ DiDi s ^ — ( x _ Xo )k! + (x + x )k 2 , 



y-ys 



DiDi. 
D s 



■<*v = (v - Vo)ki + (y + yo)k2, 



(61) 



(62) 



where (x B ,y s ) are the source positions (projected on the lens 
plane), fci and k 2 are dimensionless functions of the baryon 
mass, distances and inclinations. 

Consider the simplest case we have a face-on Kuzmin 
disk with the parameters xo = yo = 0, zq = b. For this 
axisymmetric lens, deflection points to the lens origin and 
is a function of R only. One can predict convergence and 
shear by taking appropriate derivatives. The critical (Ein- 
stein) ring can be obtained by setting x s = y s = so that 
the Jacobian A' 1 = = 0. 

Let 



dxdy 



Rm = 



4GM DiDis 
c 2 D a 



D = 



c 

ao 



6c 
Ho 



(63) 



be the Einstein ring size in GR and a characteristic MOND 
distance respectively, then TeVeS predicted critical radius R 
satisfies 



1 = 



DiDi s a(R) 



D s 



R 



2{k 3 + k 4 ), 



where 



2fc 3 = 



Ra 



R2 i + v^ 

2/c4 = (2 arcsin \) 
2k 5 ee (ttx) 



Rm DiDis 
R D s Do 
Rm DiDi s 



R DsDo 



X 



(64) 

(65) 

(66) 
(67) 
(68) 



VR 2 + b 2 ' 

where we have used ki — k 2 = ks + ki. 

The corresponding face-on Kuzmin disk plus spherical 
DM model predicts 

DiD u a(R) 



1 



D s 



R 



2(fc 3 +k B ), 



(69) 



Hence MOND can create an illusion of an isothermal halo 
lens with a terminal velocity v 2 — ^GMao through the 
terms with k^ and k&. 

Note fcs > /c4. This means we have a slightly *bigger* 
critical radius for the case with DM halo vs. the case with 
TeVeS. Likewise TeVeS would predict slightly different time 
delay between images than DM halos for a fixed Ho- These 
are somewhat surprising since the two models have identical 
rotation curves. This suggests that a combination of lensing 
and kinematics will be able to differentiate DM and TeVeS. 



8.4 MONDian effects on critical rings 

To understand the correction due to MOND, let's consider 
lensing by a spherical point lens, which can be obtained by 
letting b = 0, hence \ — L 



2fc 3 = 



Rm 
R? ' 



2fc 4 = 2fc 5 



■kRm DiDi a 



(70) 



R D D S 

ft is helpful to estimate rescaled lens-source distance in 
TeVeS 

DiDi, DiD u 



D D S 25GpcA 



< 0.1 



(71) 



for all lenses and sources. Hence 2fc 4 < 2fc 6 ~ < 
2fc 3 ~ 1 near the critical ring R ~ Rm, so the MON- 
Dian effect is never very important near critical line. This 
means critical rings (or strong lensing) can only occur in 
line of sight which pass through regimes of strong gravity 

R ~ Rm < 1 fi^- where *the MONDian effect is small*. 



9 LENSING TESTS OF TEVES 

For lenses with almost co-linear double images in the CAS- 
TLES survey, Zhao, Bacon, Taylor, Home (2006) conducted 
a detailed fit using spherical point or Hernquist profile 
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lenses. Cares have been taken in including the K-correction, 
the luminosity evolution with redshift, and the possibility of 
significant gas and extinction from dust. They applied two 
methods, using the image positions only, and using the im- 
age amplifications. They found that the mass-to-M* ratios 
calculated using the two independent methods closely agree, 
and all but two of the lenses are found to have M/M* be- 
tween 0.5 and 2. This shows that TeVeS is a sensible theory 
for doing gravitational lensing. 

However, the authors caution that there are several 
lenses in galaxy clusters which require extremely high M/L. 
Clearly more detailed models are needed, including flatten- 
ing and the cluster environment. As a first attempt in this 
direction, Angus et al. (2006) found that the lensing peaks of 
the Bullet Cluster could be explained by adding neutrinos 
in a TeVeS-like modified gravity; the phase space density 
of neutrinos at the lensing peaks requires 2eV mass to in 
order not to violate exclusion principle for fermions. In gen- 
eral gravitational lensing can be used as a useful approach 
to distinguish between theories of gravity, and to probe the 
functional form of the modification function /i. 
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10 OTHER SANITY CHECKS OF TEVES 

Sanity checks from small to large scale have also been done 
in recent papers. TeVeS is found to be 

• OK with solar system (Bekenstein & Maguijo 2006) 

• OK with Milky Way and Bulge Microlensing (no cusp 
problem, Famaey & Binney 2006) 

• Excellent description of spiral rotation curves (Mc- 
Gaugh 2005, Famaey et al. 2006) 

• OK with elliptical galaxies lenses (Zhao, Bacon, Taylor, 
Home 2006) 

• OK with galaxy clusters if with neutrinos (Angus, Shan, 
Zhao, Famaey, 2006), 

• TeVeS universe can accelerate (Zhao 2006, astro- 
ph/0610056) 

• Structures and CMB can form from linear perturba- 
tions (Dodelson & Liguori 2006). 

TeVeS is by no means a firmed established paradigm 
since many comparisons of the theory with observations are 
still unknown, but in the process of understanding and fal- 
sifying TeVeS, we learn to design clever dark matter models 
and appreciate better the robustness of GR. As it stands, 

• TeVeS is not grossly inconsistent with observations of 
lensing apart from a few outliers associated with galaxy clus- 
ters where massive neutrinos would contribute to the deflec- 
tion of the light, 

• CMB anisotropy are predictable (Skordis et al. 2005) , 

• structure formation in non-linear potential can in prin- 
ciple be followed by N-body codes (Ciotti et al. 2006) . 
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